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ABSTRACT.  The  linearized  equations  of  motion  of  a  yawing  projectile 
are  generalized  to  include  the  effects  of  wind,  and  expressed  in  a  form 
suitable  for  use  on  an  analog  computer.  They  should  be  especially 
useful  in  interpreting  the  data  obtained  with  yaw  sondes. 
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INTRODUCTION 


Reference  1  presents  an  analysis  of  the  yaw  of  certain  artillery 
shell.  The  angular  motion  had  been  measured  by  solar  yaw  sondes,  and 
was  compared  with  the  motion  predicted  by  theory.  During  the  comparison, 
the  parameters  in  the  theory,  especially  the  aerodynamic  coefficients, 
were  varied  until  good  fits  were  obtained  between  data  and  calculations. 
Similar  studies  of  what  is  essentially  the  7-caliber  Army-Navy  shaped 
projectile  are  reported  more  briefly  in  Ref.  2  and  5.  Now  digital  com¬ 
puter  programs,  in  several  versions,  are  available  which  solve  the  dynam¬ 
ical  equations  of  motion,  and  these  equations  are  complete  and  rigorous 
insofar  as  the  aerodynamic  force  system  is  known.  But  up  to  the  present 
the  turn-around  time  on  available  digital  machines  has  been  a  few  hours 
at  the  least.  Furthermore,  automatic  iterative  methods  of  parameter 
selection  appear  to  be  excluded,  for  they  require  a  starting  solution 
which  is  already  a  very  good  fit  to  the  data.  If,  for  example,  there  is 
anywhere  a  discrepancy  in  the  phase  of  oscillation  of  even  half  a  cycle, 
an  iteration  will  almost  certainly  give  no  or  a  nonsensical  result.  The 
intervention  of  human  Judgment,  guided  by  considerations  such  as  those 
steming  from  Dr.  C.  H.  Murphy's  "quasi-linear"  analysis,  has,  so  far, 
been  essential.  Hence  an  analog  computer  has  been  used,  which  computes 
the  motion  in  about  one-and-a-half  times  the  real  time  of  flight  and 
provides  concurrently  a  plot  which  can  be  directly  compared  with  the 
plotted  data.  The  limitations  of  existing  analog  devices  have  required, 
in  turn,  a  simplification  of  the  equations  of  motion. 

The  simplification  consists  in  breaking  the  problem  into  two  parts. 
The  first  part  is  the  computation  of  a  smoothed  trajectory,  essentially 
that  of  a  particle,  neglecting  oscillations.  (Spin  history  is  also 
fitted  here.)  In  the  second  part  the  transverse  angular  motion  is  com¬ 
puted  from  equations  which  have  been  simplified  to  the  extent  that  the 
spherical  trigonometry  of  the  rotation  group  representation  is  linear¬ 
ized.  On  the  other  hand,  nonlinearities  in  the  aerodynamic  force  and 
moment  system  are  specifically  included,  so  far  as  known  or  suspected. 
This  part  uses  the  velocity,  trajectory  angle,  and  altitude  (via  air 
density)  from  the  first  part.  Feedback  of  the  oscillations  into  the 
first  part  is  not  permitted.  An  apparent  exception  to  the  last  state¬ 
ment  arises  because  it  has  sometimes  been  desirable  in  the  first  part, 
to  take  account  of  the  effect  of  yaw  on  drag.  But  this  has  been  done  by 
Inserting  a  "mean  square  yaw"  as  a  smooth,  known  function  of  the  time, 
its  values  having  been  chosen  by  estimation  with  the  aid  of  a  prelimi¬ 
nary  fit  of  the  oscillations. 
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The  system  of  equations  used  so  far  makes  no  provision  for  the  ef¬ 
fects  of  wind.  As  mentioned  in  Ref.  1,  there  are  clear  indications  that 
the  system  needs  such.  In  addition  to  making  the  oscillation  fit  better 
and  more  realistically,  inclusion  of  wind  ought  to  give  a  more  realistic 
and  more  easily  fitted  particle  trajectory.  These  points  were  again 
raised  recently  during  a  visit  to  this  Station  by  people  from  the 
Ballistics  Research  Laboratory,  Aberdeen  Proving  Ground.  Now  it  was  not 
immediately  obvious  how  to  include  the  wind  terms  in  the  equations  in  a 
way  most  nearly  consistent  with  the  levels  of  approximation  used  before. 

This  note  proposes,  and  gives  a  skeleton  derivation  of,  a  set  of 
equations  which  include  wind.  Two  cautions  are  in  order:  first,  the 
equations  have  not  yet  been  tested  in  actual  application;  second,  the 
approximations  would  be  inadequate  when  yaw  becomes  large  or  the  wind 
speed  an  appreciable  fraction  of  the  projectile  speed.  Such  conditions 
could  arise,  for  example,  near  the  suianit  of  a  very  high-angle  trajectory. 


THE  EQUATIONS 


Two  comments  are  in  order.  Because  a  so-called  yaw  sonde  or  yaw 
camera  actually  measures  not  yaw  but  the  angle  between  the  projectile 
axis  and  a  line  to  the  sun,  we  must  keep  track  of  the  actual  orientation 
as  well  as  the  yaw.  This  requires,  for  the  angular  motion,  a  system  of 
equations  of  higher  order  than  that  usually  used  in  analytical  discus¬ 
sions  of  stability.  Then  we  wish  the  particle  or  reference  trajectory 
to  be  as  close  to  the  real  one  as  possible.  This  will  require  including 
wind  in  both  of  the  normal  equations.  This,  in  turn,  together  with  our 
desire  to  be  able  to  include  mean  square  yaw  effects  on  drag  means  that 
the  natural  aerodynamic  coefficients  to  use  are  K_  and  K,,  not  K_. 
and  K^,  and  furthermore,  that  the  particle  trajectory  is  to  be  calcu¬ 
lated  neglecting  the  odd  powers  of  the  yaw  as  measured  in  the  air 
mass. 

We  will  use  reduced  velocity  U  ■  V/d  and  wind  W  *  (Wind  velocity )/d 
where  d  is  reference  diameter. 

Let 

Ur  =  U-W,  u  -  |u|,  ur  «  |Ur|. 

w  *  iWj  +  J w2 

We  define  the  aerodynamic  coefficients  and,  implicitly,  the  coordi¬ 
nate  system  by  taking  P  as  a  unit  vector  along  the  projectile  axis, 
and  writing 

P  =  i  cos  y  cos  +  J  8in  y  cos  $2  +  k  8ln  ^2 

y  -  a  + 
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U  «  u(i  cob  P  cob  0g  +  J  sin  P  cob  0g  +  k  Bin  ©g) 
P  -  a  +  ©: 

0  «  f/dm  -  Jg/d 

f/d  -  -  pd5  [KDurUr  t  KLUrx(UrxP)  +  KjP(UrxP)] 

U  ■  -  Du  U  -  LU  x(U  xP)  -  Fto(U  xP)  -  Jg/d. 
r  r  r  r  r 

The  spin  angular  velocity  is  to,  and 

D  *  pd^Kjj/ra,  L  ■  pd^j/tn,  F  ■  pd^/ra. 


•  • 

Now  u  u  =  U  •  U. 

U  ■  Ur  ■  U  uf  COB  r.uur  [l  -  0(  |w|^)]  g  U  Ur 

Also  we  set  P  *  U  cos  6  /u  +  A  with  U  •  A  «  0  and  |A  I  *  sin  &  . 

r  rrr  rr  r  r 

Then  U  xP  «  U  xA  and  U  •  (U  xP)  -  U  •  (-WxA  ). 
r  r  r  '  r  r 

The  term  in  F  of  u  u  will  then  both  be  small  and  of  first  order  in 
yaw,  and  is  to  be  neglected. 

U  x(U  xP)  -  U  x(U  xA  )  -  -U  2A 
r  r  rrr  r  r 

The  term  in  L  1b,  then,  also  first  order  in  yaw  and  should  be 
neglected.  Therefore,  we  take  (neglecting  also  0^) 


u  u  „  -  Dur  u  -  g  sin  a/d 


(1) 


U  ■  u  ( i  cos  P  cos  Sg  +  J  sin  P  cob  ©2  +  k  sin  9^) 

+  u  P  cos  ©2  (-  i  sin  0  +  J  cos  P) 

+  u  Gg  (-  i  cos  P  sin  ©g  -  J  sin  P  Bin  ©2  +  k  cob  0g) 

•  •  • 

U  p  COB  ©2  “  U  ’  (-  i  Bin  p  +  J  COB  p)  »U  • 

•  • 

u  ©  ■  U  •  (-1  cob  P  sin  0g  -  J  sin  P  Bin  ©g  +  k  cos  ©g) 
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Now 

U  •  ■  U  *  €g  ■  0 

U„  •  •  -W  •  sin  P  ~  sin  a 

Ur  •  €g  «  -W  •  «2  •  v1  cos  p  sin  6g  -  wg  cos  eg  a  -  wg 

Urx(UrxP)  *  -ur2P  +  urUr  cos  8r 

[Urx(UrxP)]  •  €1  -  -ur2P  •  *x  +  ur  cos  5r  Ur  •  €x 

[Urx(UrxP)]  .  €g  -  -ur2P  •  €2  +  ur  cos  5rUr  .  €g. 

But 

P  .  €1  -  cos  02  sin  (y  -  0)  «  -  61 

P  •  €g  «  -  cos  (y  -  0)  sin  6g  cos  02  +  sin  0l2  cos  02  02  -  eg 

Or 

[Urx(UrxP)]  •  «1  jj.  \i<^1  -  e:)  +  sin  a 

[Urx(UrxP)]  •  €g  ^  -  up2(0g  -  e2)  -  urv2  . 

By  straightforward  expansion  followed  by  linearization  we  find 
(UrxP)  .  €x  ^  -u(02  -  ©2)  -  Wg 

(UrxP)  .  €g  u  u(^  -  61)  -  wx  sin  a 

Thus 

uP2  -  Du^  sin  a  +  Lur  [u  (0  -  sin  a] 

♦  Fto  [u(02  -  e2)  +  w2]  - 
but 

*  •  • 

0  •  a  +  0^  ,  cos  0  ~  cos  a  -  0^  sin  a 

We  take 

ua  *  -  g  cos  a/d  -  Du  w.  sin  a  (2) 

r  i 

Then 

«  2 

u  Lur  6lr  +  Ptuu52r  -  01  g  sin  a/d  (3») 

4 


NVIC  TP  1*898 


with 

*  $1  ~  sin  °/u  (since  u/ur  <*  l) 

52r  ■  «2  *  ®2  *  Vu 

And 

u  0 g  *  +  Lur2  62 r  “  Rou  6lr  ( 5b) 

The  Magnus  force  terns  (with  P  as  a  factor)  have  been  included  only 
for  completeness.  It  is  recommended  that  they  not  be  used. 

The  torque  is  taken  to  be 

t  «  -  pd^  K.u  cdP  +  p  d^KMu  (U  xP) 

A  r  m  r  r 

-  pd^  KTaiPx(UrxP)  -  pd^Kg^O 

where  Q  is  the  transverse  angular  velocity  and  0  •  P  i  0. 

If  I  and  I,p  are  the  polar  and  transverse  moments  of  inertia, 
we  have  ^ 

t  -A  &P“p  * 

=  1  a>  P  +  I  a>  P  +  I_  Q 
p  p  *T 

.  The  easiest  way  to  find  Q  is  to  note  that  P  *  0  xP,  and  hence 
PxP  *  Px(  Q  xP)  ■  Q  . 

We  find,  to  the  approximation  we  need 

•  • 

0  ~  02(  i  sin  y  -  J  cos  y )  +  k  y  , 

o  -  o(022  +  y2)  , 

and 

•  •  • 

P  ~  -  y  ( i  sin  y  -  J  cos y )  +  k 

But 

p  .  p  -  0 
or 

I  i  *  t  •  P 
P  " 

• 

a)  -v  -  Au  to  (k) 

-  r 
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with  A  ■  p  K./l  . 

A  p 

But 

Q  +  ffi  P  -  [t  »  P(  T  •  P)]/^  -  r'/^ 

where 

5  “  "  Ip/IT 

Neglecting  a  and  terms  In  y  ,  and  setting  *  i  sin  y  -  J  cosy  , 

••  •  , 

02-  55  0i  ^  to  at  +  T '  • 

M  • 

^1  +  5  ^2  -  ^  ‘  k^IT* 

Expanding  the  cross  products  as  before, 

»i  •  • 

<f>l  +  <5  02  *  -  Hur  +  Muru  8lf  +  T<nu  8^  (5a) 

(here  Ha  is  neglected) 

\  -  ®  ^i  *  -  Hur#2  +  Muru  62r  *  T5u  8lr 

-  ®  (fi-SSLS  +  Duf  sin  a)  (5b) 

One  should  be  able  to  drop  the  distinction  between  u  and  u  in  all 
terms  except :  r 

1.  The  drag  terms  in  Eq.  1 

2.  The  gravity  term  of  Eq.  2 
5.  Equation  4 

4.  The  main  moment  term  (that  in  M)  in  Eq.  5 

The  last  item  (4)  is  not  really  an  Inconsistency,  for  the  difference 
between  u  and  uy  will  affect  the  frequency  of  the  precessions. 

The  coordinates  of  the  particle  trajectory  are,  of  course,  obtained 
from 

x  ■  d  u  cos  a  1 

y  »  d  u  cos  a  J  (6) 

The  numbered  equations  form  the  desired  system.  The  dimensionless 
aerodynamic  coefficients  are  proportional  to  the  relative  air  density 
p/e0. 
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M,  T,  H,  L,  and  F,  if  used,  may  be  taken  as  functions  of  Mach 

2  2  2 

number  ■  U  /a,  and  of  instantaneous  6„  ■  6,  *  6^_  .  Dependence  on 

.  .  r  r  ir  2r 

^  and  $2  also  be  included  if  desired,  so  long  as  the  well-known 
symmetry  conditions  are  satisfied.  D  will,  as  mentioned  above,  depend 
on  Mach  number  and  a  smoothed  version  of  6f  ,  chosen  before  the  partic¬ 
ular  angular  motion  calculation.  A  probably  needs  to  depend  only  on 
Mach  number. 
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